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Abstract 

We study classical and quantum self-similar collapses of a massless scalar 
field in higher dimensions, and examine how the increase in the number of 
dimensions affects gravitational collapse and black hole formation. Higher 
dimensions seem to favor formation of black hole rather than other final states, 
in that the initial data space for black hole formation enlarges as dimension 
increases. On the other hand, the quantum gravity effect on the collapse 
lessens as dimension increases. We also discuss the gravitational collapse in a 
brane world with large but compact extra dimensions. 
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I. INTRODUCTION 



Since the Choptuik's discovery |IJ of critical phenomena in gravitational collapse of a 
scalar field, there have been much attention to the gravitational collapse just at the threshold 
of black hole formation. Similar phenomena to Choptuik's results were quickly found in 
other systems numerically 0- (Tl| and analytically [fL^] - |19| too,[] suggesting that they were 



limited neither to scalar field matter nor to spherical symmetry. There are two possible final 
outcomes of gravitational collapses, distinguished by whether or not a black hole is formed 
in the course of evolution. Which end state is realized depends on initial data specified by a 
control parameter p that characterizes the gravitational interaction strength in the ensuing 
evolution. For p < p* the gravitational field is too weak to form a black hole, while for 
p > p* a black hole is produced. The critical solution with p = p* just at the threshold of 
black hole formation acts as an intermediate attractor for nearby initial conditions and has 
an additional symmetry: discrete or continuous self-similarity. 

While critical gravitational collapse has been originally studied in four spacetime dimen- 



sions, there have been several attempts to study it in higher-dimensional spacetime |2T 
27| . As most attempts (including string theory) to unify gravity with other forces require 
additional spatial dimensions, it may be meaningful to consider higher dimensional collapse 
problem. Although the extra dimensions are compactified to a small size and usually unob- 
servable at low energy, the extra dimensions may be opened in the large energy limit of the 
Planck scale, which is relevant as we approach singularity in gravitational collapse. 

Moreover, in recent years there have been explosive interests in the large extra dimension 
scenario ||32| , in which the size of extra dimensions can be as big as sub-millimeter and the 
fundamental scale of the higher-dimensional gravity as low as a TeV range. An exciting 
consequence of this scenario is the possibility of production of black hole at realistic future 



colliders . The Schwarzschild radius of so produced black hole will be much less than the 
size of extra dimensions because the mass of the black hole is typically of the order of the 
higher- dimensional Planck scale of TeV. Therefore, in this scheme the gravitational collapse 
will be essentially higher-dimensional and quantum gravitational. We can also expect the 
formation of black holes with much larger mass than the fundamental mass scale, but with 
smaller Schwarzschild radius than the size of extra dimensions. Then, black holes can be well 
treated as general relativistic objects in higher dimensions. Such a black hole may be formed 
from a macroscopic initial distribution of matter prepared in four-dimensional world, but 
whether black hole are formed from the given initial data depends on the number of extra 
dimensions, because the final stage of collapse will be entirely higher-dimensional. Therefore, 
we expect that one can determine the number of large extra dimensions by investigating the 
gravitational collapse and the dimensional dependence of black hole formation, at least, 
in principle. In this sense it will be valuable to study the gravitational collapse in higher 
dimensions and the dimensional dependence of black hole formation. 

In the next section, we shall start with a brief summary of Frolov's works [fj^] , which treat 
the higher-dimensional scalar field collapse with restriction to the spherical symmetry and 
self-similarity and examine an analytic solution. We then examine the role of dimensionality 
of spacetime in the gravitational collapse and black hole formation. Higher dimensions seem 



1 See the excellent review by Gundlach [20], and more references there in 
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to favor black hole in that, as dimension increases, the critical value of the control parameter 
decreases and so the initial data space for black hole formation enlarges, while the size of 
the apparent horizon increases near the critical collapse. 

We will also consider, in Section [II, the quantum gravity effects on the collapse of the 
scalar field by extending the work of Ref. p9|j30[] to higher dimensions. Since the mass of the 
black hole near the critical collapse can be arbitrarily small, the quantum gravity effects may 
play a significant role in such collapse. We also examine the role of spacetime dimensionality 
in such quantum collapse. The quantum gravity effects on the collapse decrease as the 
spacetime dimension increases. Finally, in the last section, we discuss the gravitational 
collapse in the brane world with large extra dimensions. 



II. DIMENSIONAL DEPENDENCE OF GRAVITATIONAL COLLAPSE 

We begin this section with a short review of Frolov's solutions of Ref. [^] for gravitational 
collapse of a minimally coupled scalar field in n dimensions. The action of the system is 



S= f d n Xy/^ 
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where k = 87rM* ^ and M* is the fundamental scale of the n-dimensional gravity theory. 
The scalar field has a mass dimension of (n — 2)/2. With the spherical symmetry and the 
self-similarity ansatz 



ds 2 = e 2a{T \-x 2 dr 2 + dx 2 ) + r 2 dQ 2 n _ 2 , r = xy(r) and 
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after taking into account the surface term the Einstein-Hilbert action is reduced to 
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where the dots denote derivative with respect to r and we have introduced a dimensionless 



parameter x = x 
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where (= M t 1 ) is the fundamental length. We have integrated 



out the angular variables, and V n -2 is the volume of (n — 2) -dimensional unit sphere, i.e., 
V n -2 = / d£L n -2 = 2n( n ~ 1 ^ 2 /T[(n — l)/2]. From this reduced action, one can easily see 
that ct = constant, and after appropriate rescaling of coordinates a can be set to zero. The 
classical equations of motion for and y are then given by 
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Here, c (> 0) is a dimensionless integration constant and characterizes the initial distribution 
of the massless scalar field. It physically represents the gravitational interaction strength in 
the collapse and its inverse gives the measure of the inhomogeneity of the collapse .0 



2 The constant c does the same role with parameters A and Co of Refs. and (^8|, respectively, in 
which self-similarly collapsing null fluid and dust are examined in different spacetime dimensions. 
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From Eq. (|5]), we see that the collapsing system can be replaced by the system of a 
particle (living in y space) with zero energy under the influence of the potential 

vM = \(i-y 2 )- c ^ 2{ - 3) - (6) 

The collapse is then parallel to the motion of a particle which starts from infinity (y = oo) 
and rolls toward zero (y = 0). 

To see when a black hole is formed, we introduce the local mass function M(r, x) such 

as 

2K M(T ' X) -(Vr) 2 . (7) 



(n-2)K- 2 r»" 3 

Clearly, if n — 4, this reduces to the one usually used in four-dimensional spacetime, and if 
the spacetime is static, it gives the correct mass of Schwarzschild black holes in n-dimensions 
||3 1|| . We can then localize the apparent horizons given by (Vr) 2 = 0, which translates to 
y 2 — y 2 = 0. With this fact, from Eq.(|5|) we see that the black hole is formed if the value of 
y reaches 

yAn = c 1 ^ n - 3 \ (8) 

There are three types of solutions depending on the value of the constant c, i.e., sub- 
critical (c < c*), critical (c = c*) and supercritical (c > c*) solutions, where c* is a crit- 
ical value of c characterized by potential having second order zero, i.e., V c \ = = at 
y* = yj(n — 3)/(n — 2), and it is given by 

nn-2 
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In the supercritical case where V c \ < everywhere, the value of y reaches ?/ah and a black 
hole is formed by the collapse. On the other hand, in the subcritical case where there exists 
a region with V c \ > 0, the value of y turns around at the turning point and never reaches yAn, 
which is located in the forbidden region. In other words, the scalar field collapses, interacts 
and disperses leaving behind a flat spacetime without forming a black hole. In the critical 
case y takes infinite time to reach yxu and y = 0, and the spacetime is asymptotically flat 
but contains a null, scalar-curvature singularity. 

The local mass function of the black hole formed in the supercritical collapse is given 
from Eq.(0) by 



There, A and Co characterize the initial distributions of null fluid and dust in each collapses. Their 
inverses also give the measure of the inhomogeneity of the collapse. For large values of their inverses 
one has highly inhomogeneous collapses where the outer shells collapse much later than the central 
ones [Q. The common characters of the parameters c, A and Co can be easily observed from 
comparing the energy densities of the systems. In fact, we could choose other parameters that 
prescribe the initial data, instead of c. However, with the above physical meaning it seems to be 
more favorable to choose c than other parameters. 
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M AH = ^^rtt = (IL^, 1 ^ V5? = kz<9^W~. (10) 

Clearly, on the apparent horizon the mass becomes unbounded as x — > oo. The spacetime 
fails to be asymptotically flat for supercritical case. As a result, the total mass of the black 
hole cannot be written in power-law form in terms of initial data. This divergence is an 
unnatural artifact of self-similarity. In the self-similar solutions this kind of divergence is 
unavoidable. In order to avoid this divergence one has to consider the influx of the scalar field 
with a finite duration only. This is equivalent to introducing a cutoff x c for the x-coordinate. 
Since this is not allowed in the strict sense of a self-similar solution, we just suppose that 
any behavior of fields in the outer region beyond the boundary x c does not seriously disturb 
the self-similar evolution described by the action (Q). One may then interpret the solution 
as a nearly self-similarly collapsing (at least in the inner region of x < x c ) solution with a 
finite mass of 

M BH (n-2)c 1 / 2 T/ fx c y~ 3 



On the other hand, in Ref. fl8j a more realistic model that represents such type of collapse in 



four spacetime dimension is constructed by a cut and paste method: cutting the spacetime 
along a null hypersurface and then joining it to an exterior solution (e.g., an outgoing Vaidya 
solution). The mass of the black hole is finite and always takes the form Mbh oc (c — c*) 1 / 2 . 
We easily expect from the symmetry of the system that a simple extension of cut and paste 
method of Ref [TH] to higher dimensions leads to a similar result 

M BH cx(c-cT, (12) 

where 7 = (n — 3)/ \j2(n — 2) is the scaling exponent and it was identified in Ref. |21| . 

It is interesting to observe that the critical value c* depends only on the spacetime dimen- 
sion n, and significantly and monotonically decreases as the number of dimensions increases 
as shown in Eq.@ and Figure 1. The value t/ah a ls° increases as the number of dimensions 
increases. These facts tell that the black hole window for initial data enlarges as the space- 
time dimension increases, because the parameter c specifies the initial distribution of the 
massless scalar field. Let us consider a system collapsing from a given initial distribution 
of the scalar field specified by a value c smaller than the critical value c* in n- dimensional 
spacetime. A black hole can then be formed through the gravitational collapse from the 
initial distribution in higher than n dimensions, while with the same initial distribution 
the scalar field will disperse leaving behind just a flat spacetime without forming a black 
hole in n or lower dimensions. Thus, we can conclude that higher dimensions favor a black 
hole rather than a flat spacetime with dispersive scalar field when the massless scalar field 
self-similarly collapses. 

Even though we have considered only the self-similarly collapsing scalar field to study 
the effect of increase in dimensionality on the results of gravitational collapse, we would 
expect similar effects for the dimensional dependence of black hole formation in different 
systems from the universality of the critical phenomena. In fact, it was already observed 
in a few different systems. The authors of Ref. p6| investigated self-similar gravitational 



collapse of null dust fluid in higher-dimensional spacetimes and showed that as the number 
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FIG. 1. Values of c* and j/ah for critical collapse in various dimensions. 



of dimensions increases window for naked singularity shrinks, that is, higher dimensions 
favor black hole rather than naked singularity. This could be seen from the fact that A c 
decreases significantly as the dimension of the spacetime increases, where A c is the critical 
value of the control parameter A which specifies the initial distribution of the null fluid as 
the parameter c in our case. A c depends on the spacetime dimension n also and has the 
same value with c* as A c = (n — 3) n ~ 3 / (n — 2) n ~ 2 . [] While the final outcome of the collapse 
is a naked singularity if A < A c , it is a black hole if A > A c . This shows shrinking of 
the naked singularity window for initial data. A similar behavior has also been shown in 
the gravitational collapse of an inhomogeneous dust cloud described by higher-dimensional 
Tolman-Bondi spacetimes in Ref. |28], in that the effect of extra dimensions appears to be a 
shrinking of the naked singularity initial data space (of 4D) or an enlargement of the black 
hole initial data space, which is characterized by the parameter (o there. 



III. QUANTUM COLLAPSE IN HIGHER DIMENSIONS 

Black holes are well understood as general-relativistic objects when their mass Mbh far 
exceeds the higher-dimensional gravity scale M*. As Mbh approaches M*, the black holes 
will become quantum gravitational objects and their properties become complex. This gives 
rise to an obstacle in calculating the production of black hole or analyzing the gravitational 
collapse. Hence, in the small mass limit it would be physically interesting and important 



to study the quantum gravity effect in the collapse. From the mass formulas Eqs. (|TT|) an d 
(|I~2|), we see that there are two independent small mass limits, x c — > and c — > c*. In the 
case x c ~ £* the black hole mass is of the order M*. As c —>■ c* it vanishes, so the classical 
description at the final stage would not be valid. 

It is the purpose of this section to study quantum mechanically the black hole formation 
and to investigate how quantum gravity effects modify the classical picture and the dimen- 
sional dependence of gravitational collapse in higher dimensions. We use the Arnowitt- 
Deser-Misner (ADM) formulation to quantize the model and analyze the black hole wave 
function. In order to canonically quantize the system, we start from setting a = in the 



Here we have redefined A according to our convention, so this is different by the factor 2 ( n 2 ) 



from that appears in Ref. 26] 
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metric Eq.(|3]) and introduce a lapse function N(r) to write the metric as 

ds 2 = -x 2 N 2 (T)dT 2 + dx 2 + x 2 y 2 (T)dn 2 n _ 2 . (13) 
With this metric the action ([$]) is rewritten after taking into account the boundary term as 



S = V n . 2 £:~ 2 I d X dTNy n - ' 
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The canonical momentum densities of this system for y and are defined as 

dC n — 3 y n ~ 4 y 
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Then the action 5 can be rewritten in the ADM formulation as 



S = V r . 
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where TC is the effective superspace Hamiltonian density. The lapse function acts as a La- 
grange multiplier, so varying the action with respect to N we get the Hamiltonian constraint: 
H = 0, where 
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Following Dirac's quantization method, the Hamiltonian constraint becomes a quantum 
constraint on the wave function 



H{y,7i y ; 0,vr )|$ >= 0. 
The Hamiltonian constraint is then reduced to 
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where we have defined the canonical momenta = V n -2 

„n-2 



m-2 



Jdx^i 



Kifi, and K = 



V n -2Pl~ 2 J dx = V n -2X n ~ 2 . K represents the volume of (n — 2)-dimensional sphere with 
radius x. It has an infinite volume in the limit x — > 00, and the canonical momenta and 
the Hamiltonian are divergent. This divergence is an unnatural artifact of self-similarity as 
the black hole mass also diverges. In order to avoid this divergent Hamiltonian we have to 
consider the influx of scalar field with finite duration only, as done to avoid the divergent 
black hole mass in the previous section, and we introduce a cutoff x c in the x-coordinate. 
This cutoff is related to a scale of the inner region where the assumption of self-similarity 
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is valid. Since the large x c limit is the large black hole limit, it will be equivalent to the 
classical limit. 

When we quantize the system, the first term of equation ([19]) contains an ordering 
ambiguity between operators y and tl y . For many choices of the factor ordering, the effect of 
the factor ordering can be parameterized by a constant a, and the corresponding Hamiltonian 
is obtained by the substitution of U 2 — > — y~ a [d y y a d y ). However, in this paper we do not 
adhere to the generality of the operator ordering, and for simplicity we take the simplest 
ordering, i.e., a = 0. Of cause, due to this-like free fixing of the operator ordering the 
validity of our result may be limited. However, we expect that the factor ordering does not 
affect any of the semiclassical calculations done in this paper.0 We will let the choice of a 
be dictated by convenience. The Wheeler-De Witt (WDW) equation ([H]) then becomes 



'd_Y + (n-2)(n 
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where K c = V n -2^c 2 ■ The wave function can be separated into the scalar field and gravi- 
tational field parts as 



V(<j ) ,y) = e ±in ^<5>(y)- (21) 
Then the Wheeler-De Witt equation reduces to one-dimensional Schrodinger type equation 



2 [dy, 
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where the quantum mechanical potential VwDwiy) is given by 
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with c related to the conjugate momentum density tt^ by c = [{n — 2)/(n — 3)}kti 2 . As easily 
observed from Eqs. (f|) and ([15|) c is the same parameter with that used in the previous 
section. This potential is different from its classical counter part (^) by the multiplication 
factor [(n — 3) 2 K 2 / ' K 2 ]y 2<yn ~ 4 \ The shape of the potential in each dimension depends on two 
parameters, the initial data c and the cutoff x c . As in the classical collapse, c determines 
the supercritical (c > c*), the critical (c = c*) and the subcritical (0 < c < c*) collapses. 
The cutoff x c makes the top of the potential shift away from zero. 

In the classical collapse, the collapsing system was replaced by the system of a particle 
living in y space with zero energy under the influence of the potential V c i, and the collapse 



4 For arbitrary a, taking dy = y a d y the WDW equation ( |22| ) is written as [— (l/2)dy + 
Vwdw]®(X) = 0, where Vwdw(X) = V 2a VwDw{y)- On the other hand, as we will see soon after 
in the subsection 3. A, the black hole formation rate depends on the factor a = J dy \/\V\yDw(y)\ ■ 
If we would calculate it with the semiclassical solution &(Y) instead of &(y), then the rate would 
depend on a = J dY ^J\VwDwiX)\- However, we can easily see that a = a. Therefore, the factor 
ordering does not affect any of the semiclassical calculations. 
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parallels to the motion of a particle which starts from the infinity (y = oo) and rolls toward 
zero (y = 0). On the other hand, the Wheeler-De Witt equation ( P2f) looks like a Schrodinger 
equation, and wave functions which satisfy the WDW equations describe a particle living in 
y space under the influence of the potential Vwdw- Therefore, from experiences in quantum 
mechanics it seems to be plausible to replace the quantum collapse with the scattering 
problem of a wave function living in y space. The boundary condition for the wave functions 
for black hole formation then is such that the wave function should be incident from the 
spatial infinity in y space, i.e., y = oo, and some part of them be reflected by the potential 
barrier back to the infinity of y but the remaining part be transmitted toward the black hole 
singularity inside apparent horizon. With armed these boundary conditions, we compute 
the rate for black hole formation. 

Unfortunately, it is impossible to find exact solutions of the Wheeler-De Witt equation 
( p2f) except for n = 4. Therefore, we will try to see some aspects of quantum collapse by 
using the WKB method in several limiting cases: i) x c 3> and c C c*, ii) c ~ c*. Here, 
x c ^> means that x c is large enough for the top of the potential to be far from zero, so 
that the WKB is valid near the top of the potential. 



A. Subcritical collapse: x c ^> and c ^ c* 



In the case of subcritical collapse, the potential Vwdw m any dimension has two zeros which 
we will call y n \ and y n 2 (> y n i)- We then define three separated regions, I (0 < y < y n i), II 
{yni < y < yra) and III (y > y n 2)- In this regime we assume c is small enough compared to 
c* and x c large enough, so that the semi-classical analysis valid. Since we are interested in 
the collapse of the scalar field starting from y = oo, we assume that the wave function has a 
purely incoming flux toward the singularity at y = in region I, and the current of incident 
wave function in region III is unity for convenience. With these boundary conditions, the 
WKB wave functions in each region are given by 
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Except for the phase factor, the wave function fl24j) depends only on a which determines the 
tunnelling and reflection amplitudes of the waves. The wave function $ni indeed satisfies the 
boundary condition for black hole formation, in that the first term describes the incoming 
component and the second term describes the outgoing (reflected) component at the far 
region (y > 1). 



S 



The existence of nonvanishing wave <3>i in region I implies that there is a possibility to 
form a black hole quantum mechanically, even if it is classically forbidden. We now calculate 
the rate of a black hole formation in the subcritical collapse. The incoming flux at spatial 
infinity computed with Eq . (j26|) is 



j m (y) = Im <$>* m {y)in y <S> ixi {y) 



1. 



(28) 



Some part of the incoming flux is reflected by the potential barrier back to the infinity 
y = oo but the remaining part transmits toward the black hole singularity at y = 0. On the 
other hand, the reflection rate at the infinity is similarly found as 



Pr = — = tanh 2 a. 

Jin 



(29) 



The rate for black hole formation is the ratio of the transmitted flux to the incident flux. 
From the flux conservation we obtain the rate for black hole formation 



r~T — — 

Jin 



I — — — sech 2 a. 



(30) 



obtained by a direct computation 



This transmission rate coincides with the rate near y 
using the transmitted wave Eq . (pl|) . 

As expected previously, the classical limit is obtained in the large limit of K c (or equiv- 
alently x c ). In a spacetime of dimension n, a depends only on the cut off as cx (x c /£*) n ~ 2 . 
As a increases the transmission rate Pp exponentially decreases and vanishes, while the 
reflection rate Pr approaches to unity. Therefore, most part of the scalar field bounces back 
without forming a black hole in the limit of large x c . 
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FIG. 2. The rate for black hole formation in diverse dimensions for various values of the 
control parameter c and a value of the cutoff x c = 1.6£pi. 



On the other hand, the value of a increases as c decreases, and has a maximum value 
when c = for a given cutoff x c and dimension n. The rate for black hole formation falls 
off as c decreases, and closes to the minimum rate in the limit c — > 0: 
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~t min — >• seen ; , as c — * 0. 31 

* 2r (f ) 

Since the parameter c characterizes the gravitational interaction strength in the collapse 
and small c means small interaction strength, it is reasonable that the rate for black hole 
formation falls off as c becomes small. However, it seems surprising that there exists a mini- 
mum rate for black hole formation with vanishingly small gravitational interaction strength 
of limit c — > 0. This is not an artifact of our approximation by WKB method. We can also 
observe this aspect in Ref. ||30|| , in which the transmission rate is exactly calculated using the 
analytic solution of Eq. (p2[). Unless c = definitely, the potential is singular at the origin 
at y = 0, and the incident waves are not reflected at the origin.^ Note that the curvature 
singularity at the origin acts like a sink for the incident waves. Therefore, the transmission 
rate does not become zero in the limit c — > 0. The minimum rate may have a different form 
from Eq. (|3T1), because in the limit c — > the potential becomes very flat near y = y n \ 
and is very close to zero in a broad range. So our WKB approximation would not be valid 
near y = y n \. The numerical works in Fig. § confirm the nonvanishing transmission rate 
in the limit c — > 0. On the other hand, we have to consider the position of the apparent 
horizon. If the apparent horizon is at y = as in the classical collapse, then there is no black 
hole formation despite of the existence of the transmission flux in the region I. However, 
the apparent horizon and black hole mass will be affected by quantum effect. The quantum 
mechanical effect will shift the position of the apparent horizon from y = to y ~ 1/ K% as 



discussed in Ref. |3Q 



The rate for black hole formation doesn't seem to vary simply with spacetime dimensions 
for a fixed initial data c and cutoff x c . The n dependence of a is complicated as can be 
seen from Eq.(^7]). The volume of (n — 2)-dimensional unit sphere V n -2 included in K c 
does not monotonically behave but has a maximum value at n — 7, while the other factors 
monotonically increase as n increases. However, with some generic values of the cutoff x c , a 
behaves simply and Pt monotonically decreases within four to eleven dimensions, in which 
we are interested from the unified theory perspective like string theory. Since the integral in 
Eq.(^7|) is not explicitly calculable in general dimensions, we present numerical values of Pt 
for several different c in Figure 2. The result shows that the rate for black hole formation 
by quantum tunneling falls off as the number of dimensions increases for a given initial data 
c and fixed cutoff x c . 

B. Near-critical collapse: c ~ c* 

In this regime we cannot directly apply the usual WKB approximation, because the height 
of the potential is too low and the derivative of the potential is almost zero near turning 
points. So the semiclassical solution such as $n of previous case is not valid near the top of 



5 The case c = should be distinguished from the limit c — > 0. If c = definitely, the potential 
is not singular but is zero at y = 0. The incident waves are totally reflected at the origin, so the 
transmission rate is then zero. Notice that the spacetime described by the metric Eq. (|2|) does not 
have a singularity, but is regular at the origin. 
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the potential. Note that WKB is valid only when the condition SX/X <C 1 is satisfied, where 
A is the wave length of the wave functions. The spirit of the WKB method is that one find 
a solution of the wave equation near the turning point by approximating the potential with 
a linear function. Then one matches the solution with the other two semiclassical solutions 
outside the transition region. 

We develop an approximation which has a common spirit with the usual WKB method. 
We approximate the potential Vwdw with a quadratic function near the position of the 
top and find an analytic solution of the wave equation with the quadratic potential. We 
then match the solution with the other two semiclassical solutions in regions I and III 
choosing appropriate coefficients. This approximation method is not so different from the 
conventional WKB method but only a slight deformation from the WKB method, besides 
that we introduce a quadratic approximation instead of a linear approximation. So it has 
much common sprit with the usual WKB method, and it will give approximately good 
answer as much as the usual WKB method does. 

We represent the near critical limit by setting c = c*(l + 5) with \S\ <C 1, where the 
collapse is supercritical if 5 > and subcritical if 5 < 0. The Wheeler-De Witt potential is 
then separated as 



(n - 3) 2 K 2 
2n 2 y 2 



V N c{y) = Krit(y) - 



-c*S, 



(32) 



where V crit (y) is potential for critical collapse and can be written as 



V CI it(y) = 



(n - 3) 2 K t 



2k 2 



( 2 *2\ ; 

[y -y ) 



n-A 



k + 1 y *2(n-k~5) y 2k-2 



k=0 



n-2 



(33) 



Krit(y) has a second order zero at y* = y (n — 3)/(n — 2) as the classical potential V c \ does. 
Near y = y*, the potential may be approximated as 



Vnc 



(n - 3fK 2 



2n 2 y 2 



c*5 + \4* 2{n - 5) (y 2 - y* 2 ) 2 



2 V 2y 2 ' 



(34) 



where 



V- 

v+ 

E 



in ~ Z?K 2 C 2(n _ 5) 



{n - S) 3 K 2 2(n _ 3) 

.2 y 



2 k 2 
{n - 3) 3 K 2 



y 



2k 2 



1 + 2C * 6 „*-2(n-3) 



n 



(35) 

(36) 
(37) 



In fact, there could be several other choices in choosing the approximate quadratic potential. 
But we have chosen above approximate potential upon the following advantages: First, with 
this choice the WDW equation is exactly solvable. Second, it well approximates the full 
potential near the origin at y = as well as around the top, so the solutions for the 
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approximate potential behaves like $/ in the region I. Lastly, in 4- dimensional case, it 
corresponds to the exact potential, not approximate one. 

With above approximate potential, the Wheeler-De Witt equation allows solutions in 



terms of the confluent hypergeometric function [35) : 



$ NC (y) = De~ l V^y 2 / 2 (iJv^y 2 Y M(a,b,iJv^y 



where D is a constant and 



- 1 l n 
/i =4~2 g ' 

b = 1-iQ, 



(38) 



(39) 



Elf), 



with Q = (Ey* 2 v 2 — 1/4) 1 / 2 , v = V~/(Ey* 2 ). There exists another independent solution of 
the Wheeler-De Witt equation with the potential (j34]) . However, we will not write down 
that solution, because the solution ( |3"%| ) is enough to describe the formation of a black hole. 
Remind that we follow the analogy to the scattering problem of quantum mechanics and 
the boundary conditions mentioned previously. According to the boundary conditions, it 
seems reasonable to choose only purely ingoing wave functions in region I. The same wave 
function was used to describe the four-dimensional formation of a black hole by a scalar field 
collapse |I0 . 



We now match this solution with the semiclassical solutions in regions I and III. In fact, 
the solution ( |38"D is also valid near y = 0, because the approximate potential Vnc has the 
same asymptotic form as the full potential Vwdw- Therefore, we only need to match the 



solution fl4"T| ) with the semiclassical potential in region III. Using the asymptotic form of the 
confluent hypergeometric function 



lim M(a,b,z) « r(6) 

|z|^oo 



r(6 



+ 



e z z a-b 

T(a) 



7T 37T 

- - < arg(2;) < — , 



(40) 



we get the asymptotic form of the wave function at large y: 

T(b) 



i\JV + y 2 f YV^ 2 /2 + 



T(a) 



i\/V+V 2 ) e^M, (41) 



where a — /i — 1/4 — iy/Ey* /2 and b — a — /i = 1/4 + i\/Ey* /2. The first term describes 
the incoming wave function and the second term describes the outgoing wave function. 
Comparing the asymptotic form (f4"T|) with the semiclassical potential in region III of the 
form 



J J y p(y')dy' 



+ 



D-2 



/p(y) \Jp(y) 

we obtain the asymptotic form of the wave function in the region III 



(42) 



V m {y) w D 



- a—fi 



p 2 (y) 



-if y P (y>W , £W 

r(a) 



- b—a—fi 



p 2 (y) 



J J y p(y')dy' 



(43) 
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where p = \/2Vwdw- Note that the asymptotic solution ( |4~ID is of the same form as the 
semiclassical solution (|4"2]) with p = ^—2Vnc- 

Then the rate for black hole formation is easily obtained from the two solutions (IO) and 



p *E ^ f( ^ +Q) S inhvrg 

Jin COShf(^*-Q) 

It is easy to check that this transmission rate coincides with that obtained by using current 
conservation at the spatial infinity 



p = 3jc = co&%(Q + y/Ey*) c _ nQ = 1 _3* (A5) 

R Jin COshf(Q- y/Ey*) " Jin' 

This determines the probability for black hole formation in the near critical collapse for 
values of |5| <C 1. 

In the large cutoff limit x c ^> it will be easier to see the behavior of the rate Pp. In 
the large x c limit the rate for black hole formation should be calculated separately for 5 > 
and 5 < 0. For 5 > 0, Eq.(H) becomes asymptotically 



exp 



y/26K c /n-3\i" 



\Jn — 3k \n — 2, 

In the case 5 < 0, the asymptotic value of Eq. (fP]) in the large x c limit is 

V2SK C fn-3\ n/2 ~ 

~ exp 



(46) 



\fn — 3k \n — 2 



(47) 



This is the tunneling rate for black hole formation in the subcritical collapse near the critical 
case. The case that 5 = has to be considered separately. The asymptotic form is then 

e - 5 i^VZ+^rzyl) sinh 77/ _ 1/4 

-. (48) 



coshf(V^+ y/V- - 1/4) 

As in the previous case, the classical limit is obtained in the very large cutoff limit x c ^> £*. 
The rate Pt exponentially depends on K c . In the very large x c limit Pt becomes unity for 
5 > 0, which implies that a black hole is formed from the collapsing scalar field as in the 
classical case. For S < 0, Pj becomes zero. This tells that the scalar field bounces back 
without forming a black hole. For 8 = 0, Pp becomes 1/2. 

On the other hand, the dimensional dependence of the rate Pp is more complicated, 
because the volume factor V n ^2 included in K c does not behave simply but reaches to the 
maximum value at n = 7, while the product of all the other factors in the exponent of Eq . ([46|) 
monotonically increases with fixed c and x c as dimension n increases. However, for generic 
values of x c the absolute value of the exponent monotonically increases between four to 
eleven dimensions. As the spacetime dimension n increases the rate for black hole formation 
Pt increases for 5 > 0, while Pr decreases for 5 < 0. Therefore, the rate for black hole 
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-0.4 -0.2 0.2 0.4 

FIG. 3. The rates for black hole formations in near critical collapses in diverse dimensions for 
a value of the cutoff x c = 1.5£p[. 

formation by quantum tunnelling falls off as dimension increases in the subcritical collapse 
with 5 < 0, while black hole formation is easier in higher dimensions for 5 > 0, as shown 
in Fig. [5| This implies that the aspect of quantum collapse becomes close to the classical 
aspect as dimension increases. This can be interpreted as the quantum gravity effects on 
the result of gravitational collapse lessens with increase of the spacetime dimensionality. 



IV. DISCUSSIONS: GRAVITATIONAL COLLAPSE IN A BRANE WORLD 
WITH LARGE BUT COMPACT EXTRA DIMENSIONS 

We have considered self-similarly collapsing system of a massless scalar field in higher 
dimensions, both classically and quantum mechanically. We examined the dimensional de- 
pendence of such gravitational collapse and formation of black hole. Higher dimensional 
spacetime seems to favor formation of black hole rather than other final outcomes, in that 
the initial data space for black hole formation enlarges as dimension increases, i.e., the crit- 
ical value c* of the control parameter which describes initial distribution of the scalar field 
falls off as dimension increases. 

We have studied black hole formation quantum mechanically by the collapsing scalar 
field. We canonically quantized the system and found the wave functions by using the WKB 
method. The results of quantum collapse obtained in four- dimensional spacetime go 



over to higher- dimensional quantum collapse. The wave functions for black hole formation 
have not only incoming flux but also outgoing flux in both the subcritical and supercritical 
cases. The classical evolution of black hole formation is recovered in the limit of large cutoff 
x c . On the other hand, the rate for black hole formation by quantum tunnelling in the 
subcritical collapse lessens as spacetime dimension increases, that is, it becomes harder to 
form a black hole in higher dimension in the subcritical collapse. The reflection rate in the 
supercritical (but close to the critical) collapse, however, decreases as dimension increases: 
The black hole formation is easier in higher dimensions in the supercritical collapse. This 
implies that the quantum gravity effect on the gravitational collapse lessens with the increase 
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of spacetime dimension. 

We now discuss the gravitational collapse in a brane world with large but compact extra 



dimensions first explored by Arkani-Hamed et. al. p2 |. The dimensional dependence of the 
final result in gravitational collapse implies that if we are living in the brane world with large 
extra dimensions, then the aspect of gravitational collapse is drastically different from that 
expected in usual four- dimensional world. Even though we prepare a macroscopic initial 
distribution of matter in four-dimensional spacetime, which final product is resulted from 
the collapse will depend on the number of extra dimensions, because if the Schwarzschild 
radius is smaller than the size of extra dimensions, the collapse will be higher-dimensional 
rather than four- dimensional at the final stage of the collapse. 

Suppose that one prepares an initial distribution of the scalar field, of which typical 
size is much larger than the extra dimension size, so that the collapse is effectively four- 
dimensional at the early stage of the collapse. Here, we assume that the three-brane has 
negligible self-gravity of its own as in usual large extra dimension scenarios and so the brane 
can be identified as a surface of vanishing extrinsic curvature. Therefore, if we restrict to 
black holes of size tbh much less than the size of extra dimensions L, then the geometry near 



the black hole will be very well approximated by a n-dimensional Schwarzschild solution [31 



If a black hole is formed from matters on the brane, the symmetry requires that the brane 
pass through the equator of the black hole. If the initial distribution is prepared (taking 
c > c 4 and large enough x c ) to form a black hole with Schwarzschild radius tbh larger than 
the extra dimension size L (rBH 3> L), then the collapse will be effectively four- dimensional 
event to an observer outside the horizon, because the observer can not see the collapse after 
the apparent horizon forms. The geometry of the resulting black hole is simply a product of 
four- dimensional Schwarzschild solution and a (n — 4)-dimensional torus, and its dynamical 
properties are determined by four- dimensional physics. 

However, if the initial distribution is prepared (if a black hole is formed) so that the 
resulting black hole have Schwarzschild radius of <C tbh "C L, the collapse goes under 
the size of the extra dimensions and becomes essentially n-dimensional at the final stage of 
the collapse.^ Therefore, the final result of the collapse can be different from that expected 
in the purely four- dimensional world. If one prepares the initial distribution specified by 
c > C4, then a black hole is formed as expected in four- dimensional world because c 4 > c* . 
Since, for r 3> L, the geometry of the resulting black hole is well approximated by the 
four- dimensional Schwarzschild black hole and the mass measured on the brane is the same 
as the mass in the bulk, it may look like a usual four- dimensional Schwarzschild black hole. 
However, it has different thermodynamic properties from those of a purely four- dimensional 



6 The mass of the black hole can still be much larger than the four-dimensional Planck scale 
and can be treated as a general relativistic object. In the usual large extra dimension scenario 
[32], the size of extra dimensions L is related to the effective four-dimensional Planck scale Mp\ 
and the n-dimensional gravity scale M* by L ~ M~ 1 (Mp\/M if ) 2 ^ n ~ 4l \ while the Schwarzschild 
radius of n-dimensional black hole is tbh = M~ 1 \1Q-kM-qy{/ (n — 2)V n -2M^] 1 ^ n ~ 3 \ So the black 
hole mass of size r B H ~ L is M BH ~ M*(Af P i/M*) 2 (™- 3 )/(™- 4 ). If M* is of order a TeV, then 
M BH ~ io 32 ("- 3 )/("- 4 )M* ~ 10 32(n-3)/(n-4)-2i g xhis is much larger than the four-dimensional 
Planck mass, e.g., for n = 6, M BH ~ 10 27 g and, for n = 10, M BH ~ 10 16 g, about the mass of the 
Earth and the typical mass of primordial black holes, respectively. 
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black hole ||33|| . The Schwarzschild radius is larger than it would be for a four- dimensional 
black hole. This means that the temperature is lower, the evaporation rate is slower, and the 
horizon area (the entropy of the black hole) is larger. If c < c* , then there is no black hole 
formation and the scalar field disperses back to the four-dimensional world. The spectrum of 
the dispersed scalar field will not be like that in a purely four-dimensional collapse, because 
the very light Kaluza-Klein modes with masses starting at ~ 1/L couples to the matter and 
some fraction of the energy radiate off of the brane into the bulk. 

The case that the initial distribution is prepared to have the control parameter of c* < 
c < c\ may be more interesting. Then, an n-dimensional black hole will be formed and there 
will be no dispersing scalar field. This is not the one expected if we are living in purely four- 
dimensional world. If our world is purely four-dimensional, or the extra dimension is very 
small and of the order of quantum gravity scale as taken in usual higher-dimensional unified 
theories, then there will be no black hole formation and the scalar field will disperse back 
to the four-dimensional world. Therefore, we can expect that one can probe the existence 
and determine the number of large extra dimensions by examining the gravitational collapse 
and the dimensional dependence of black hole formation, (at least) in principle. 

Of course, we have over-simplified the problem. At the final stage when the collapse ef- 
fectively occurs in n dimensions, the collapse is not like the n-dimensional spherical collapse 
discussed in the previous sections. The collapse is not spherically symmetric and self-similar, 
even if we have prepared the initial distribution to be spherically symmetric and the col- 
lapse to occur self-similarly at the beginning, at which the collapse can be described by the 
four-dimensional general relativity. The collapse may look like that of a disk-like rather 
than a spherical distribution at the late stage, because the matter does live only on the 
three-brane. Therefore, the discussions in the previous two paragraphs suffer from lack of 
physical meaning. Despite of this poor physical meaning, the purpose we presented the dis- 
cussions upon the over-simplified assumption is simply to illustrate the probing mechanism 
of the large extra dimensions. In this sense, the discussions in the two previous paragraphs 
seems to have its necessity. For rigorous treat, we have to consider the collapse of disk-like 
distribution of matters and the critical phenomena. However, it seems impossible to find 
analytic solutions for collapses of disk-like matter distributions. Some numerical simulations 
for such collapses may be helpful. But, this is beyond the scope of this paper. 

In spite of all, we easily expect that the brane world scenario makes interesting predictions 
about the formation of small black holes. In particular, the existence of the large extra 
dimensions would affect on the formation of primordial black holes, which are formed as a 
result of the density fluctuation in the early universe. If there was large extra dimensions, 
then mass spectrum of the primordial black hole would not be like that predicted in purely 
four-dimensional theory. It will be interesting to investigate their detailed phenomenology. 
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